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ABSTRACT 

We present a systematic analysis of the rotation curves of 187 galaxies with masses 
greater than 10 ^°Mq, with atomic gas masses from the GALEX Arecibo Sloan Sur¬ 
vey (GASS), and with follow-up long-slit spectroscopy from the MMT. Our analysis 
focuses on stellar rotation curves derived by fitting stellar template spectra to the 
galaxy spectra binned along the slit. In this way, we are able to obtain accurate ro¬ 
tation velocity measurements for a factor of 2 more galaxies than possible with the 
Ha line. Galaxies with high atomic gas mass fractions are the most dark-matter dom¬ 
inated galaxies in our sample and have dark matter halo density profiles that are 
well fit by Navarro, Frenk & White profiles with an average concentration parameter 
of 10. The inner slopes and of the rotation curves correlate more strongly with stel¬ 
lar population age than with galaxy mass or structural parameters. At fixed stellar 
mass, the rotation curves of more actively star-forming galaxies have steeper inner 
slopes than less actively star-forming galaxies. The ratio between the galaxy specific 
angular momentum and the total specific angular momentum of its dark matter halo, 
Rj, correlates strongly with galaxy mass, structure and gas content. Low mass, disk- 
dominated galaxies with atomic gas mass fractions greater than 20% have median 
values of Rj of around 1, but massive, bulge-dominated galaxies have Rj = 0.2 — 0.3. 
We argue that these trends can be understood in a picture where gas inflows triggered 
by disk instabilities lead to the formation of passive, bulge-dominated galaxies with 
low specific angular momentum. 

Key words: galaxies:haloes; galaxies: formation; galaxies: structure ; galaxies: stellar 
content 


1 INTRODUCTION 

Rotation curves are an important probe of the distribution 
of both the luminous and the dark matter mass in galaxies. 
They have been studied since the 1970’s primarily as a way 
to constrain the density profiles of the dark matter halos that 
are thought to surround all galaxies. Because dark matter 
increasingly dominates the mass distribution at large dis¬ 
tances from the centre of a galaxy, most observational stud¬ 
ies have focused on rotation curves measured using the HI 
line at 21 cm wavelengths. This is because the atomic gas in 
galaxies generally extends significantly further out than the 
stars. Many galaxies, however, have holes in the HI distribu¬ 
tion at small radii, so if one wishes to measure the shape of 
the dark matter potential over a wide range of scales, the HI 
data needs to be complemented by rotation measurements 
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from optical Ha or CO mm line data (see Sofue & Rubin 
2001; Combes 2002 for comprehensive reviews). 

There is a significantly smaller body of work on the 
statistical properties of rotation curves. Persic et al. (1996) 
compiled 1100 rotation curves of spiral galaxies and claimed 
that a single functional form with galaxy luminosity as 
the only parameter could provide a “universal” description. 
The most luminous galaxies show slightly declining rotation 
curves in the outer part, following a maximum in the disk. 
Intermediate galaxies have nearly flat rotation curves across 
the disk. The least luminous galaxies have monotonically 
increasing rotation curves across their whole disk. These 
results led Persic et al (1996) to conclude that the dark- 
to-luminous mass ratio increases with decreasing galaxy lu¬ 
minosity and mass. Later studies of early-type disc galax¬ 
ies (SO-Sab) by Sofue et al (1999) and Noordermeer et al 
(2006) revealed large discrepancies with the rotation curves 
predicted by the Persic et al. fitting formula. These authors 
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found that the shape of the rotation curve depends on the 
bulge-to-disk ratio of the galaxy; galaxies with concentrated 
light distributions and larger bulges reach their maximum 
rotation velocities at smaller radii than galaxies with smaller 
bulges and more diffuse light distributions. 

One limitation of past statistical studies is that because 
they rely on Ha, HI or CO line data, they are biased to 
galaxies with high gas content and/or ongoing star forma¬ 
tion. In principle, galaxy rotation curves can also be derived 
from stellar absorption lines, but this requires high signal- 
to-noise spectra and accurate fitting of stellar templates in 
order to accurately measure shifts in velocity along the slit. 
In recent years, this has become possible thanks to the in¬ 
creasing efficiency of optical spectrographs and the avail¬ 
ability of high resolution stellar template spectra spanning 
a wide range in wavelength (e.g. Martinsson et al 2013a). 

Moran et al. (2012) presented long-slit spectroscopy of 
174 galaxies drawn from the GALEX Arecibo Sloan Digi¬ 
tal Sky Survey (SDSS)(GASS; Catinella et al. 2010), which 
was designed to measure the neutral hydrogen content of a 
representative sample of ~ 1000 galaxies with stellar masses 
larger than lO^^M© and redshifts 0.02 < 2 < 0.05 uniformly 
selected from the SDSS (York et al. 2000) and Galaxy Evolu¬ 
tion Explorer (GALEX; Martin et al 2005) imaging surveys. 
GASS observations detect HI down to a gas fraction limit 
of ~ 3%. A companion project on the IRAM 30m telescope, 
GOLD GASS (Saintonge et al 2011), has obtained GO (1-0) 
line measurements and molecular gas masses for a subset 
of these galaxies. In addition, a variety of stellar popula¬ 
tion and structural parameters derived from GALEX/SDSS 
imaging and spectroscopy are available for the full sample. 

This paper presents a systematic analysis of the rota¬ 
tion curves derived from the final sample of 236 longslit 
spectra obtained over the period from October 2008 to April 
2012 using both the Blue Ghannel Spectrograph on the 6.5m 
MMT telescope and the Dual Imaging Spectrograph on the 
3.5m telescope at Apache Point Observatory (APO). As part 
of the procedure for measuring SFR and gas-phase metal- 
licity profiles for the galaxies in the sample, Moran et al 
(2012) derived both stellar rotation curves from fitting tem¬ 
plate stellar absorption line spectra to the galaxy spectra 
binned along the slit, and Ha rotation curves from fitting 
Gaussian profiles to the Ha emission line after subtraction 
of the best-fit spectrum. 

As we will show, the agreement between the Ha and 
stellar rotation measurements are generally good to within 
10-20% in galaxies where both can be accurately deter¬ 
mined. We are able to obtain accurate stellar rotation veloc¬ 
ity measurements for a factor of 2 more galaxies than is pos¬ 
sible with Ha alone (187 out of 236 galaxies), permitting a 
much wider exploration of galaxy parameter space. The stel¬ 
lar rotation curves generally do not extend beyond ~2.5Re 
and thus do not in general sample the flat part of the ro¬ 
tation curve with high accuracy. Nevertheless, we will show 
that we can still use them as probes of, a) the inner mass 
distribution of the galaxy, b) the radius at which dark mat¬ 
ter begins to dominate over baryons (stars and gas), c) the 
angular momentum content of the galaxy as parametrized 
by the spin parameter A. After presenting our data analy¬ 
sis methodology in Section 2, we examine how these three 
quantities depend on a variety of galaxy properties in Sec¬ 
tion 3. In Section 4 , we select a subset of galaxies where dark 


matter dominates well within the optical radius, and we ex¬ 
amine constraints on the shape of the dark matter density 
profile. In section 5, we summarize and conclude. We adopt 
a Hubble constant of Ho = 70 km s“^ Mpc”^ throughout 
the paper. 

2 DATA ANALYSIS METHODOLOGY 

2.1 Rotation curve estimation and sample 
selection 

The procedure for measuring rotation curves is described in 
detail in Moran et al (2010). Working outwards from the 
galaxy centre, each spectrum is binned spatially to achieve 
a signal-to-noise of at least 6 (A“^) in each bin, the mini¬ 
mum needed for a reliable velocity estimate. For each binned 
spectrum, the effective spatial position is determined by cal¬ 
culating the luminosity-weighted average radius of all spatial 
positions that entered the co-add. 

The radial velocity as a function of radius is measured 
in two ways: 1) by cross-correlating each spectrum against 
superpositions of template spectra from Bruzual & Ghar- 
lot (2003) stellar population synthesis models to determine 
the velocity, 2) by fitting directly for the centroids of the 
Ha emission lines. The methods yield generally consistent 
curves at radii where we can fit templates and measure emis¬ 
sion line centroids. This is shown in Figure 1, where we show 
plots of the difference between the velocity determined by 
fitting stellar templates (Vabs) and the velocity infered from 
the shift of the centroid of the Ha line {Vna)- The velocity 
differences are normalized by dividing by 0.5(1468 + Vna) to 
yield a fractional difference. We compute the average frac¬ 
tional distance for each galaxy with where Vabs and Vhc are 
both well-measured over at least 10 spectral bins, and we 
plot this quantity as a function of different global galaxy 
parameters in each panel of Figure 1. 

As can be seen, Vhc is generally larger than 146s, but 
this is a small effect (10% systematic effect on average). A 
similar effect was found by Martinsson et al (2013a), which 
these authors attributed to “asymmetric drift”, i.e. the fact 
that stars can be supported against gravitational collapse 
by both rotation and by random motions. If this is the case, 
one might expect the average velocity difference to correlate 
strongly with the bulge-to-disk ratio of the galaxy, and for 
the velocity difference to get systematically larger in the in¬ 
ner regions of the galaxy. We can explore this hypothesis by 
plotting the velocity difference as a function of S* for each 
spectral bin, where E* is the local stellar surface mass den¬ 
sity of stars measured from our derived stellar mass profiles. 
This is shown in Figure 2. Black points show the fractional 
difference between Vna and Vabs for individual spectral bins, 
while the red squares show the average value of the fractional 
offset in fixed bins of E*. As can be seen, the average frac¬ 
tional offset increases from 0.03 at stellar surface densities 
less than 10 ®Mq kpc“^ to 0.15 at stellar surface densities 
of 10® Mq kpc“®, which is consistent with the “asymmet¬ 
ric drift” hypothesis. The bottom panels of Figure 1 show 
that the largest velocity diferences averaged over the rota¬ 
tion curves of individual galaxies occur for galaxies with red 
NUV-r colours and high inclinations. This might be caused 
by dust extinction effects or systematic errors in the mea- 
surment of the centroid of the Ha when the signal-to-noise 
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Figure 1. The mean fractional difference between the veloc¬ 
ity estimated by fitting stellar templates to the longslit spectra 
iYabs) t)y fitting to the centroid of the Ho emission line {Vho.) 
is plotted as a function of different global galaxy parameters in 
each panel: stellar mass (top left), r-band concentration parame¬ 
ter (top right), global NUV-r colour (bottom left), and axial ratio 
(bottom right). 



(V.bs-V„J/ 0.5(V^,^+V„J 


Figure 2. The fractional difference between the velocity esti¬ 
mated by fitting stellar templates to the longslit spectra ( 1 ^( 163 ) 
and by fitting to the centroid of the Ha emission line (Vhoi) is 
plotted as a function of local stellar mass surface density S* Black 
points show results for individual spectral bins, while red squares 
show the average value of the fractional difference as a function 
of S*. 


is low. In either case we expect the stellar rotation measure¬ 
ments to be more reliable. In what follows, we will analyze 
the stellar rotation curves without attempting to correct for 
asymmetric drift or extinction effects, which would require 
more detailed modelling. 

We proceed to inspect each of the 236 stellar rotation 
curves by eye and we select those that have measured veloc¬ 
ities with fractional errors less than 20% across a significant 
portion of the optical disk and those that are regular enough 
to fit with a smooth function of the form 

V{R)^Vma.R/{R°‘+r':y^°‘ + AV ( 1 ) 

(Bohm et al. 2004; Moran et al 2007,2010), where R is the 
radius, a and Va are free parameters that govern the shape 
of the rotation curve, and AV is the offset of the galaxy’s 
central velocity from the redshift obtained from the SDSS 
spectrum (also left as a free parameter). This results in a 
sample of 187 galaxies, of which 89 (i.e. around a half) also 
have Ha rotation curves measured over the same radii as 
the stellar rotation curves. The properties of this sample are 
shown in Figure 3. The black histograms show the stellar 
mass, concentration index and NUV-r colour distributions 
of the selected sample of 187 galaxies, while the red his¬ 
tograms are for the galaxies where the rotation curves were 
either very poorly measured or too irregular to work with. 
As can be seen, the selected sample of 187 galaxies with 
good stellar rotation curves spans a wide range in stellar 
mass, concentration and colour; only the very most mas¬ 
sive (M* > IO^^Mq), reddest (NUV-r> 6) and most con¬ 
centrated (i? 9 o/i? 5 o > 3.25) galaxies are systematically ex¬ 
cluded. 



Figure 3. The black histograms show the distribution of the 
number of galaxies as a function of stellar mass (left), concentra¬ 
tion parameter (middle) and NUV-r colour (right) for our selected 
sample of 187 galaxies with well-measured, regular stellar rota¬ 
tion curves. The red histograms show the distributions for the 
49 galaxies where stellar rotation curves could not be accurately 
determined, or where the rotation curves were too irregular to fit 
a function of the form given in the text. 


2.2 Parameters derived from the rotation curves 

Our stellar rotation curves are well-determined out to radii 
of between 7 and 20 kpc. This is generally not far out enough 
in the disk to measure the shape of the outer rotation curve, 
but we are able to accurately measure the inner slope and 
determine whether the gravitational potential is dominated 
by baryons (stars-fgas) or by dark matter over the optical 
radius of the disk. In addition, as we will discuss, we measure 
the total angular momentum content of the stellar disk, as 
parameterized by the spin parameter A. 
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(i) The inner slope is defined to be the slope of the best- 
fit linear relation between Vabs and R, evaluated from the 
center of the galaxy out to a radius of ±rs, where Vs is the 
fitting parameter in equation (1). The inner slope provides 
a measure of the shear in the central regions of the galaxy. 
In order to be able to compare massive galaxies with less 
massive galaxies in a meaningful way, we scale the veloci¬ 
ties by dividing by Vmax, where Vmax is an estimate of the 
maximum rotation velocity of the galaxy. We also scale the 
radii R the radii by dividing by Rso , where R50 is the radius 
containing half the total r-band light of the galaxy. Because 
the errors on the fitted value of Vmax in equation ( 1 ) are 
large for galaxies where the outer regions of the rotation 
curve are not very well sampled, we estimate Vmax using 
the baryonic Tully-Fisher relation given in Catinella et al 
( 2012 ) (logVmax = 0.2371oga;-0.251;a; = M* -b IAMhi). 

(ii) The dark-matter domination radius Rso{DM) is de¬ 
fined as the radius where the enclosed dark matter mass 
exceeds the total baryonic mass. The total mass enclosed 
within radius R can be estimated directly from the rotation 
curve after correction for inclination effects. The inclination 
of the galaxy is estimated using its r-band axial ratio b/a 
as cosi = h/a. The total baryonic mass enclosed within ra¬ 
dius R is the sum of the stellar and gas mass within this 
radius. Stellar mass-to-light profiles for each of the galaxies 
in our sample are derived by fitting SDSS 5-band optical 
photometry to stellar population synthesis model grids us¬ 
ing standard techniques (see for example Salim et al 2005). 
The mass-to-light ratio will depend on the assumed initial 
mass function (IMF). For simplicity, we have assumed a 
fixed Chabrier (2003) form. 

Because we only have single-dish measurements of the 
21cm and CO (1-0) lines, the gas mass within radius R must 
be approximated using assumed profiles. It has been known 
for some time that there is a tight relation between the char¬ 
acteristic size of the HI disk D1 and the total HI mass in 
the disk (Broeils & Rhee 1997), 

log Mhi = 1.96 log Dl + 6.52, ( 2 ) 


where Dl is defined as the diameter where the face-on cor¬ 
rected angular-averaged HI column density reaches 1 Mq 
pc“^ (corresponding to 1.25 x 10^° atoms cm“^). This rela¬ 
tion has recently been demonstrated to extend up to total 
HI masses of 2 x IO^^Mq and disk diameters of 100 kpc 
(Wang et al 2013). Wang et al (2014) show that a simple 
analytic expression of the form 


T.Hi{r) 


h exp(-r/rs) 

1 -b /2 * exp(—r/rc) ’ 


( 3 ) 


where I\, I 2 , Vg and are free parameters, can provide 
an excellent description of the HI density profiles of nearby 
gas-rich and normal spiral galaxies. Bigiel & Blitz (2012) 
show that total (atomic-bmolecular) gas profiles are well- 
characterized by a function of the form 

T.gas/T.traasit = 2. le" (4) 


where J^transu is the radius where the atomic gas surface 
density equals the molecular gas surface density, and R 25 is 
the optical radius of the galaxy. 

We construct atomic and molecular gas surface density 
profiles for the galaxies in our sample as follows. Equation 
(2) allows us to specify Dl, given the observed total HI mass 


of the galaxy. Figure 6 of Wang et al (2014) shows that 
rs/Rl = rs/(0.5Dl) lies in the range 0.15 to 0.3. Also based 
on this data, we allow r^ to range from O.OSrg to O.Trg. Fi¬ 
nally, we also require that f 27rEHi(r)rdr = Mhi- We first 
generate a family of HI profiles that satisfy these constraints. 
If the galaxy has detected molecular gas, then we add in 
molecular gas at each radius such that equation (4) is satis¬ 
fied and we only keep those profiles where the integral over 
the resulting molecular gas profile yields a value close to the 
observed total H 2 mass. Finally we add each acceptable gas 
profile to the stellar profile and compute Rso{DM) as the 
value averaged over each S*(r) -b combination. In 

practice, because stars dominate the baryonic mass budget 
for the majority of the galaxies in our sample, the uncer¬ 
tainty in Rsq{DM) arising from the variety of acceptable 
gas profiles is small, typically 10-20 %. 

(iii) The spin parameter A. Observational measurements 
of the specific angular momentum in galaxies place con¬ 
straints on the way galaxies acquire their mass and an¬ 
gular momentum. The specific angular momentum is often 
parametrized in terms of a dimensionless spin parameter A. 
Here, we adopt the definition of galaxy spin parameter given 
in equation (2) of Dutton & Van den Bosch (2012) 


\gal — 


{Jgal/Mgal) 
V2R vir K ir 


( 5 ) 


where Jgai is the total angular momentum of the galaxy, 
Mgai is the total baryonic mass of the galaxy, and R^ir and 
Vvir are the radius and circular velocity of its dark matter 
halo. As discussed by Dutton & Van den Bosch (2012), if 
Xgai can be measured for a representative set of galaxies, the 
angular momentum ratio Rj, defined as the ratio between 
the galaxy specific angular momentum and the specific an¬ 
gular momentum of the halo. 


Rg = 


{Jgal/Mgal) 

{Ja ir /Ma ir') 


^gal 

Xhalo 


( 6 ) 


can then be constrained. The distribution of Xhaio has been 
extensively studied using N-body simulations (e.g. Bett et 
al 2007) and semi-analytic models of disk galaxy formation 
(e.g. Fu et al 2010, Guo et al 2011) typically assume that 
Rj — 1 , because this leads to good agreement with the ob¬ 
served sizes of galaxies. In this paper, we examine the valid¬ 
ity of this assumption. 

Calculating the total angular momentum for each galaxy 
in our sample, given its rotation curve and baryonic mass 
density profile is straightforward: 


Jgai^ 211 j {T,t,{r)-\-T,gas{r))Vrot{r)rdr, (7) 

Accurate estimates of R^ir and Vvir cannot be made for indi¬ 
vidual galaxies in our sample. Instead, we adopt the parame¬ 
terized relation between galaxy stellar mass and dark matter 
halo mass at z = 0 derived using the galaxy abundance¬ 
matching technique in Moster et al (2013), 

M^Mhaio = 2N [{Mnaio/Mi)-^ + (Mhaio/Mi)^] , ( 8 ) 

with N = 0.035, Mi = 3 x lO^^M©, j} = 1.5 and 7 = 0.6. 
The virial radius Rvir can then be computed from the re¬ 
lation Mhaio = 200pcrit(4n'/3)i?®ir, where pcrit is the crit¬ 
ical density of the universe, and K,ir is given by Vvir = 
{GMhalo/Rvir)^^^- 
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3 RESULTS 

3.1 Inner slope and core radius trends 

In Figure 4, we show plots of the inner slope (with veloc¬ 
ity scaled by Vmax as described in the previous section, and 
radius scaled by the half-light radius of the galaxy R 50 ), as 
functions of a variety of different global galaxy parameters. 
Black points show inner slope values measured for individ¬ 
ual galaxies, while red squares show the running median. 
We have binned the points along the x-axis such that there 
are always 12 galaxies in each bin, and the errorbar on the 
median value of the slope is calculated by boot-strap resam¬ 
pling the galaxies in each bin. 

As can be seen from the first two panels, there is no sig¬ 
nificant trend in inner slope as a function of the stellar mass 
or stellar velocity dispersion of the galaxy. The inner slope 
also does not appear to vary as a function of galaxy struc¬ 
tural parameters such as the stellar surface mass density 
or the concentration index of the light. The inner slope is, 
however, quite strongly correlated with parameters describ¬ 
ing the age of the stellar populations and the gas content of 
the galaxy. We note that the detailed behaviour of this trend 
is different for each of the 5 parameters we plot. A roughly 
continuous trend is seen as a function of the global NUV-r 
colour, with the bluest galaxies exhibiting inner slopes that 
are twice as steep as the reddest galaxies. For other parame¬ 
ters such as atomic and molecular gas mass fraction, specific 
star formation rate and 4000 A break strength, no depen¬ 
dence is seen for the gas-poor/old galaxy population, but a 
dependence does set in at higher gas mass fractions and at 
young ages. The 4000 A break index D„(4000) and the spe¬ 
cific star formation rate SFR/M, are both measured within 
the 3 arcsecond SDSS fibre aperture and hence probe the 
stellar populations in the very central regions of the galaxy. 
As can be seen in the last two panels, there is again a de¬ 
crease in inner slope for galaxies with the very youngest 
central stellar populations. 

In Figure 5, we plot our galaxies in the two-dimensional 
planes of stellar surface mass density, concentration index, 
NUV— r colour and 4000 A break strength versus stellar 
mass, and we colour-code each point according to the value 
of the inner slope: cyan points correspond to inner slope 
values of less than 0.05, black points to inner slopes in the 
range 0.05 to 0.2, and magenta points to inner slopes greater 
than 0.2. In agreement with the results shown in Figure 4, 
the galaxies with small and large inner slopes separate most 
clearly in NUV-r and D„(4000). However, in this plot we 
also see that at fixed colour, stellar surface density or con¬ 
centration index, galaxies with steep inner slopes are more 
massive than galaxies with flat inner slopes. This plot helps 
us to understand results in the literature (e.g. Persic et al 
1996), which strongly emphasize the dependence of rotation 
curve shape on galaxy luminosity or mass. The strong trend 
with colour will only be found in samples like ours that in¬ 
clude weakly star-forming and early-type galaxies. 

A quantity that is closely related to the inner slope is 
the radius Ts, defined in equation (1). This is the radius 
where the rotation curve flattens away from the inner power- 
law inner region. We note that our definition of the inner 
slope of the rotation curve is mainly controlled by varia¬ 
tions in the quantity Vs/Rso, because AVc/Vmax is gener¬ 
ally ~ 1. Ta can be regarded as the “core radius” of the sum 






log M. (MJ 


log M. (MJ 


Figure 5. Galaxies are plotted in the two-dimensional planes of 
stellar surface mass density, concentration index, NUV—r colour 
and 4000 A break strength versus stellar mass. Each point is 
colour-coded according to the value of the inner slope: cyan points 
correspond to inner slope values of less than 0.05, black points to 
inner slopes in the range 0.05 to 0.2, and magenta points to inner 
slopes greater than 0.2. 


of the dark matter halo and stellar density profiles. Our re¬ 
sults imply that blue, star-forming galaxies have total mass 
distributions with smaller core radii than red, passive galax¬ 
ies. We note that previous work (e.g. Kormendy & Freeman 
2004; Spano et al 2008; Donato et al 2009) has focused on 
the inferred core radii dark matter halos after subtracting 
the contribution from the stars and gas. As we have already 
discussed, most of the galaxies in our sample are gravita¬ 
tionally dominated by baryons in their inner regions, so we 
will not attempt such a decomposition for the whole sample, 
only for the subset which are dark matter dominated within 
Rso. 


3.2 The dark-matter domination radius 

We are able to follow rotation curve out to large enough 
radius for 103 out of the 187 galaxies in our sample to mea¬ 
sure the radius Rsq{DM) where the dark matter first con¬ 
tributes more than half the mass interior to that radius. 
The properties of this subset of galaxies are shown in Fig¬ 
ure 6 . Black histograms show distributions of stellar mass, 
stellar surface mass density, concentration index, 4000 A 
break strength, colour and atomic gas mass fraction for all 
the galaxies where we are able to estimate Rso{DM), while 
red histograms show the distributions of these properties for 
galaxies where we do not trace the rotation curve far enough 
to make an estimate of this quantity. As can be seen, the low¬ 
est mass galaxies with M* < 2 x 10 ^° Mq are missing from 
the sample with measured Rso{DM), but otherwise the sam¬ 
ple covers a wide range in colour and structural properties. 

Figure 7 shows how a variety of galaxy parameters cor- 
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Figure 4. The inner slope of the rotation curve (see text) is plotted as a function of (from left-to-right and from top-to-bottom) a)stellar 
mass M*, b) stellar velocity dispersion measured from the SDSS spectrum, c) stellar mass surface density /r,, d)concentration index, e) 
NUV-r colour, f)molecular gas mass fraction, g)atomic gas mass fraction, h)4000 A break strength, i)specific star formation measured 
within the SDSS fibre. Black points show inner slope values measured for individual galaxies, while red squares show the running median. 
We have binned the points along the x-axis such that there are always 12 galaxies in each bin, and the errorbar on the median value of 
the slope in calculated by boot-strap resampling the galaxies in each bin. 


relate with Rso{DM). Once again, we find weak correlations 
with galaxy mass and structural parameters and stronger 
correlations with galaxy colour and gas content. The two 
bottom panels of the figure show the parameters that corre¬ 
late best with the dark matter domination radius: the global 
NUV-r colour and the atomic gas mass fraction. As can be 
seen, the majority of galaxies with atomic gas mass fractions 
greater than 0.1 and NUV-r colours less than 4 become dark 
matter-dominated within 2 R 50 . Galaxies with less than a 
few percent atomic gas and NUV-r> 5 become dark matter 
dominated at radii larger than 2 R 50 . 


3.3 Spin parameter 



In Figure 8, we show plots of the differential (left) and cu¬ 
mulative (right) distributions of derived spin parameters for 
all 187 galaxies with well-measured rotation curves in our 
sample. We are able to use the full sample, because in gen¬ 
eral, the stellar rotation curves can be accurately determined 
over most of the optical radius of the disk. As can be seen, 
the median value of A is 0.01 and the Tier range in A ex- 


Figure 6. Black histograms show the distribution of properties 
of the subset of 103 galaxies for which we are able to estimate 
the dark matter domination radius, while red histograms show 
the distribution of the properties of galaxies where the rotation 
curve does not extend far enough. 
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Figure 7. The stellar masses, concentration indices, NUV-r 
colours and atomic gas mass fractions are plotted as a function of 
the dark matter-domination radius R^q(DM). Black points show 
parameter values for individual galaxies, while red squares show 
the running median of each parameter at fixed R^qi{DM). We 
have binned the points along the x-axis such that there are al¬ 
ways 12 galaxies in each bin, and the errorbar on the median is 
calculated by boot-strap resampling the galaxies in each bin. 


tends from 0.003 to 0.03. The observed median is a factor 
3-4 lower than the median value of the spin parameters of 
dark matter halo measured in N-body simulations of struc¬ 
ture formation in a ACDM cosmology (e.g. Maccio et al 
2007) and the \a scatter is about twice as large as found 
in simulations. Our results are consistent with the conclu¬ 
sions presented in Dutton & Van den Bosch (2012). These 
authors inferred spin parameter distributions indirectly us¬ 
ing ensembles of bulge-|-disk-|-halo models that were tuned 
to fit observational constraints such as the Tully-Fisher re¬ 
lation, the disk size versus stellar mass relation, and the 
relation between stellar and dark matter halo mass. Based 
on this analysis, Dutton & Van den Bosch (2012) inferred 
that the average value of A for disk galaxies was 0.019 and 
that galaxy spin did not correlate with halo mass. The big 
advantage of deriving the spin parameter directly from the 
observed rotation curve of the galaxy, is that we are able to 
examine correlations of A with different galaxy parameters. 

In Figure 9, we plot A as a function of galaxy stellar 
mass, concentration index, atomic gas mass fraction and 
4000 A break index strength. Once again, black dots show 
results for individual galaxies, while the red line shows the 
running median. The spin parameter correlates with all of 
the plotted parameters. Low mass, low concentration, gas- 
rich galaxies with young stellar populations have larger spin 
parameters than high mass, high concentration, gas-poor 
galaxies with old stellar populations. This is not unexpected 
since it is believed that the formation of bulge-dominated 
(early-type) galaxies involves processes such as interactions 
or mergers, which cause gas to lose angular momentum (e.g. 




log X log X 

Figure 8 . The differential (left) and cumulative (right) distri¬ 
butions of the spin parameters A derived for the 187 galaxies in 
our sample with well-measured rotation curves. 




Figure 9. The spin parameter A is plotted as a function of stellar 
mass, concentration index, atomic gas mass fraction and 4000 
A break strength. Black points show spin parameter values for 
individual galaxies, while the red line shows the running median 
of the spin as a function of these parameters. 


Barnes & Hernquist 1996) Interestingly, galaxies with the 
lowest masses and highest atomic gas mass fractions have 
a median spin parameter of ~ 0.03 — 0.04, which is much 
more consistent with the predictions of disk formation mod¬ 
els where the gas conserves its angular momentum. In addi¬ 
tion, we infer from Figure 9 that the scatter in spin param¬ 
eter for galaxies that are similar in stellar mass, structural 
parameters and stellar population properties is considerably 
smaller than the scatter in A for the population as a whole. 
Larger samples are needed for accurate estimates of the scat¬ 
ter in A for galaxies with fixed properties. 


4 DARK MATTER HALO PROFILE FITS 

In two seminal papers, Navarro, Frenk & White (1996,1997) 
introduced the concept of a “universal” profile (hereafter 
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NFW profile) that was able to describe the dark matter 
density profiles of halos with masses ranging from those of 
dwarf galaxy halos to those of rich galaxy clusters, inde¬ 
pendent of the initial density fluctuation spectrum and the 
values of the cosmological parameters. The NFW profile has 
the form 

pK = (9) 

Pent (r/r,)(l + r/r,)2’ 

where Vs = r 2 oolc is a characteristic radius and pertt = 
3H^/8 ttG is the critical density {H is the current value of 
the Hubble constant); 5^ and c are two dimensionless param¬ 
eters, commonly termed the characteristic overdensity and 
concentration. The virial mass and radius of the halo are de¬ 
fined through the relation M^ir = 200pcT.it(47r/3)i?®ij., and 
5e and c are then linked by the requirement that the mean 
density within is 200pcT-it. 

Rotation curve fitting is one important way to test 
whether the NFW profile can fit the dark matter halos of 
real galaxies. There are many papers in the literature with 
seemingly contradictory results, but the most recent consen¬ 
sus appears to be that the NFW profile is a reasonably good 
fit for luminous galaxies (De Blok et al 2008; Martinsson et 
al 2013), but is a poorer fit for low surface brightness dwarf 
systems, where the dark matter strongly dominates over the 
baryonic mass (e.g. Swaters et al 2003). One problem with 
the interpretation of the results for dwarf galaxies is that 
because the potential well depth is shallow, feedback effects 
from bursts of star formation during the formation of the 
dwarf can in principle affect the dark matter density profile 
(e.g. Pontzen & Governato 2012). 

In this section, we perform NFW halo fits to a subset of 
12 galaxies in our sample where the dark matter domination 
radius Rc,o{DM) < O.TRso- In the Milky Way, dark matter 
dominates over the baryons at radii larger than ~ 10 kpc, i.e. 
at around 3 R 50 . The galaxies that we analyse are thus much 
more dark-matter dominated than our own Galaxy and as 
shown in Figure 7, they also tend to have young stellar pop¬ 
ulations and to be more gas-rich. A detailed analysis of one 
such a galaxy in our sample, UGC8802 (GASS 35981), was 
presented in Moran et al (2010). In spite of its high atomic 
gas content (2.1 x 10^°M©), the star formation surface den¬ 
sity in this galaxy is low {T,sfr = 0.003Mq yr~^ kpc“^) 
and spread evenly across the galaxy. A sharp drop in metal- 
licity in the outer disk was found in the galaxy, which was 
later shown to be a common feature of many of the Hl-rich 
galaxies in the sample (Moran et al 2012). 

In Figures 10, 11, and 12, we present the NFW halo 
fits for our 12 selected dark matter-dominated galaxies. The 
measured values of the the rotation velocity are plotted as 
black circles. The dotted red curves show the contribution 
from the observed baryons (stars-l-gas), while the solid red 
curves show the best fit rotation curve derived from adding 
the contribution from a NFW halo to the baryons. Q The 
parameters of the best-fit NFW model are indicated in the 
bottom-right corner of each panel. 

As can be seen, the rotation curves of these 12 galaxies 
are regular for the most part, and the NFW halo-|-observed 

^ We have not attempted to correct for “adiabatic contraction”, 
the gravitational effect of the central baryonic contribution on the 
dark matter. 





R (kpc) 
GASS 19918 



Figure 10. NFW halo fits to the rotation curves of four dark- 
matter dominated galaxies. The measured values of the rotation 
velocity are plotted as black circles. The dotted red curves show 
the contribution from the observed baryons (stars+gas), while 
the solid red curves show the best fit rotation curve derived from 
adding the contribution from a NFW halo to the baryons. The 
Hl-profile weighted value of Vrot measured from the Arecibo spec¬ 
trum is plotted as a black star in each panel, while the value 
inferred from the rotation curve is plotted as a red star. 


baryon models provide reasonably good fits to all the galax¬ 
ies in our sample. There are no cases where the model fails 
catastrophically and in general the inner power-law region 
of the rotation curve is well-reproduced. We note that un¬ 
like many other studies, we allow no freedom in the as¬ 
sumed mass-to-light ratios for the baryonic component; as 
described in Section 2, these are derived from fits to the 
u,g,r,i,z spectral energy distributions, assuming a univer¬ 
sal Chabrier (2003) initial mass function. As an additional 
check, we compare the rotational velocity derived from the 
width of the 21cm line (see Catinella et al 2012 for details) 
with the velocity width predicted from weighting the rota¬ 
tion curve by our assumed HI profiles (see section 2). The 
value measured from the Arecibo spectrum is plotted as a 
black star in each panel, while the value inferred from the 
rotation curve is plotted as a red star. As can be seen, the 
two estimates generally agree to within 10-20% and confirm 
that our estimates of Vvtr are reasonably accurate. All the 
galaxies in our sample are found to reside in quite mas¬ 
sive dark matter halos: the smallest (GASS 3591) has a 
virial velocity of 185 km/s, corresponding to a virial mass 
of 2 X IO'^^Mq. The average value of the halo concentration 
index c is found to be 9.6 and the 10th-90th percentile range 
is from 6 to 14. This is in excellent agreement with predic¬ 
tions from ACDM cosmological simulations (see Figure 3 of 
Maccio et al (2007)). 

In summary, we conclude that NFW halos provide 
a good fit to the inner halos of the most dark matter- 
dominated galaxies in our sample. 
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Figure 11. Same as Figure 10, but for 4 more dark-matter dom¬ 
inated galaxies. 


GASS 51899 


GASS 4048 





Figure 12. Same as Figure 11, but for 4 more dark-matter dom¬ 
inated galaxies. 


5 SUMMARY AND DISCUSSION 

The main empirical results of our rotation curve study are 
the following: 

• The inner slopes of the rotation curves correlate more 
strongly with measures of stellar population age than with 
galaxy mass or structural parameters. At fixed stellar mass, 
more actively star-forming galaxies have steeper inner slopes 
than less actively star-forming galaxies. Another way to 


phrase this result is that the ratio of the radius Vs, where 
the rotation curve transitions away a power-law, to the half- 
light radius of the galaxy, is smaller in star-forming galaxies 
than in passive galaxies. 

• Galaxies with higher atomic gas mass fractions tend to 
be more dark-matter dominated than galaxies with low HI 
content. 

• The distribution of spin parameters of the galaxies in 
our sample shows that the average value of the spin param¬ 
eter of the galaxies in our sample is ~ 0.01, which implies 
that the angular momentum ratio Rj, defined as the ra¬ 
tio between the galaxy specific angular momentum and the 
total specific angular momentum of its dark matter halo, 
is less than 1. However, Rj is found to correlate strongly 
with galaxy mass, structure and gas content. The lowest 
mass, disk-dominated galaxies with atomic gas mass frac¬ 
tions greater than 20% have median values of Rj of around 1. 
The highest mass, gas-poor, bulge-dominated galaxies have 
median values of Rj of 0.1-0.2, i.e. such galaxies have sub¬ 
stantially lower specific angular momenta than their host 
halos. 

• The dark matter halo density profiles of the most dark- 
matter dominated galaxies in our sample are well fit by 
NFW models. The average concentration parameter is ~ 10, 
in good agreement with the predictions of N-body simula¬ 
tions of structure formation in a ACDM Universe. 

One might first ask whether all of these trends are qual¬ 
itatively consistent with a picture in which galaxies form by 
gas cooling and condensation within dark matter halos in 
a canonical ACDM cosmology. In such a picture, the rea¬ 
son why Rj would be less than 1 in some galaxies could 
be because gas had lost angular momentum and had been 
channelled inwards in these systems. This inward flow of gas 
would result in a build-up of higher density gas in the inner 
regions of the disk, which would then rapidly cool to low 
temperatures and form bulge stars. Gas that cools and ac¬ 
cretes from the surrounding halo is used up more quickly if 
gas inflows in the disk are common. Frequent inflows would 
thus be expected to lead to a galaxy that is passive and 
bulge-dominated. This is a possible explanation of our find¬ 
ing that the spin parameter is smaller in bulge-dominated 
galaxies with low gas content and old stellar populations. 

Gas inflows can arise in two ways: 1) when the disk is 
perturbed by a merger or interaction with another galaxy, 
2) when the disk becomes unstable and forms a bar. The 
Toomre (1964) disk stability parameter Q is larger in galax¬ 
ies with higher epicyclic frequencies and larger shear (i.e. 
larger dVrot/dR). We find that the inner slope of the rota¬ 
tion curve is large in galaxies that are still forming stars, 
which is consistent with the idea that star formation can be 
sustained for longer periods if the disk is more stable. 

The next step is to test these ideas quantitatively by 
comparing our results with N-body -|- hydrodynamical sim¬ 
ulations of disk galaxy formation in a ACDM Universe in 
which the cooling of gas in dark matter halos and its subse¬ 
quent dynamical evolution in disks is treated in a physically 
self-consistent framework. Well-resolved, dynamically con¬ 
sistent simulations rather than simplified analytic models of 
disk formation in a hierarchical Universe are also required 
to predict the effects of feedback processes from supernova 
and accreting black holes on the inner mass distributions of 
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galaxies of different types and star formation rates. They 
will also be very helpfnl in developing methodology for cor¬ 
recting for asymmetric drift. 

There are also hints of more complex dependences of 
rotation cnrve shape on star formation history. In Figure 4, 
we found that galaxies with the very youngest central stel¬ 
lar populations have rotation curves with shallower inner 
slopes. This hints that supernovae feedback processes oc¬ 
curring during strong central starbursts may be able to alter 
the inner mass distribution of galaxies. It has been proposed 
that such mechanisms may be key for resolving the observed 
rotation curves to dwarf galaxies with the dark matter halo 
profile predictions of the Cold Dark Matter model. Dynam¬ 
ical studies of larger galaxy samples observed using integral 
field spectrographs spanning a wider range of masses and 
redshifts will shed more light this issue. 
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